9 The Harmonic Oscillator 



We know that L 2 (M) is Hilbert space and it is separable, so it MUST have a complete 
orthonormal basis. For L 2 ([— it, it]) we look that the 2ir periodic solutions to the equation 



dx 2 



/ = A/ 



and these generated such a basis. Now, for L 2 (IR) we have the Harmonic Oscillator. This 

is the operator 



H 



d 2 



dx 2 



we want to solutions to Hu = Xu: 



dx 2 



+ x 2 )u = Xu, ue S(R) 



To do this, consider the two operators "Annihilation" and "Creation" 

d 



A 



dx 



x 



C 



d 

— — + x 

dx 



Now, a computation to establish at least one eigenfunction 



d 2 



dx r ' 



+ x 2 ) e-* 2 / 2 



= —x e 



dx 



so e x2 l 2 is an eigenfunction for A = 1. 

Now look at the operation of the operators A and C on each other 



ACu 



d \ ( d 
~r + x ) [ -— + x | u 

dx j \ dx 



d 2 u du 



du 



+ X- h M — X- h X U 



dx 2 dx ' dx 
So we have the relation (when we include the computation for CAu) 



d 2 



dx 2 



+ x 2 + 1 \u 



CAu=(H-l)u, ACu={H+l)u 



so 

[A, C]u = (AC - CA)u = 2u 
we already know that He~ x ^ 2 = e~ x2 ^ 2 , and Ae~ x ^ 2 = and 

HCu = (CA + l)Cu = CAC + Cu = C{AC + l)u = C{H + 2)u 
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So we now know that if u e S(R) then 

HCu = C(H + 2)u 

Set u = e~ x2 / 2 then we get the following relations by induction: 

(H - l)e- x2 ' 2 = 
(H -3)(Ce x2/2 ) =0 

H(C j e- x2/2 ) = (2j + l)C J e- x2/2 

Proposition. T/ie functions 

&e- x2 ' 2 



\Cie-* 2 / 2 \ 

are a complete orthonormal basis of L 2 ' 



First 0e x2 ^ 2 G L 2 (IR) because 0e x2 ^ 2 = Pj(x)e x2 l 2 ^ p^ (unnormalized) Hermite 
Polynomials. And by induction 



C j+1 e~ x2/2 = C(C j e~ x2/2 ) = C( Pj e- x2/2 ) = (x Pj - p'^ 2 ' 2 + (2xp 3 + p',)e 



■x 2 /2 



So pj = 2 J 'x- 7 '+ lower order terms. 

NB po, . . . ,Pn span all polynomials of degree < N, because for a polynomial qw, 



c 



In = ^Pn + Qn-i 

and so on (just pick first term then pick second...) 
Now we need to show orthonormality, that is 

( ( Pj,<Pk) = / PjWkdx = / (pj(p k dx = 0, j^k 

[H — (2j + l)](pj = 0, so we write ipj = ^j^H^j, and the above becomes 

L Vkdx = L ("^ + x2lp ) Vkdx 

these are all Schwartz functions, so we cut off at infinity and integrate by parts, 

n^oo J_ N dx ax 

r N d 

= Jim / ~ 4><P 

iV->oc J _ N dx 
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(basically, this shifts different at ion to the left) 
then we get 



<Pj -~r^ + x Vkdx = / tpjHtpkdx 



But H(p k = (2k + l)(p k so the above is just 

2k + 1 



2j + l 

then 



ifijipkdx 



(<Pn,¥k) = , , -i (<Pj,<Pk) => {<Pj,<Pk) = 



2k + 1 
2j + l 

Now we have to prove completeness, i.e. in u G L 2 (M) and J R u(fj = 0, Vj, then -u = 0. 



In fact if we have shown this then 



/ w J Cj</?j J (irr = and / up N e x ^ 2 dx = 

for all polynomials, because we know that combinations of the (fj are all polynomials. Can 
you show that this implies u = (a.e. in L 2 (IR))? Suppose that u has compact support (so 
its non-zero in a bounded set) then 

/N 
(ue~ x l 2 )p(x)dx = 0, Vpolynomials 
-N 

But the polynomials are dense in L 2 . (Think about it, this implies completeness of our 
problem) 

So any / G L 2 (R) can be written as the sum of its "Fourier Series" 



oo „ 



IPe- 2 /2| 



There are many consequences of this. For example if / G L 2 (IR) the solution to the equation 
Hu — f is 

i=o J + Jr 
In some weak sense Hu = f. Its "weak" because we don't really k now about differentiation. 
So our weak formulation is really that 

/ u (~ + x2 ^) dx = j fv> Vv? G S ( R ) 
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Completeness of solutions of the harmonic oscillator. 3ip , ipi, ip 2 , . . . in S(M) with ip the 
ground state, which means that Hip = ip , so ip = ^ exp(— x 2 /2), \\<po\\2 = 1, <fj we define 

as 



\\^ 3t Po\\L 



and 



span{(p , ...,<p n } = {p(x)e x2/2 \p poly deg(p) < N} 
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